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1 Introduction 



The inverse scattering transform method (ISM) discovered in 1967 has proved 
to be a powerful tool to construct exact solutions and to solve the Cauchy 
problem for a large variety of nonlinear integrable models of mathematical 
physics. But real physical applications are usually related to mathematical 
models with boundary conditions. For this reason the problem of adopt- 
ing the ISM to boundary value problem as well as to initial boundary value 
(mixed) problem is very important. During the last two decades this field 
of research has been intensively studied. It becomes clear that only special 
kinds of boundary conditions preserve the integrability property of the equa- 
tion given. Different approaches were worked out to look for such classes 
of boundary conditions based on Hamiltonian structures pQ , on higher sym- 
metries [2J, [3J, [3], and the Lax representation [5], [6]. Integrable initial 
boundary value problems on a half-line (in 1 + 1 case) or a half-plane (in 1+2 
case) for soliton equations nowadays is a rather studied subject. Analytical 
aspects have been developed in [TJ, [8], [9], [10], [11] where large classes of 
solutions were constructed. However, boundary value problem for the ellip- 
tic soliton equations or initial boundary value problem for regions with more 
complicated boundary are still much less investigated (see, [12] , [32] , [33] )• 

If the boundary conditions are not consistent with the integrability prop- 
erty of the equation then the standard version of the inverse scattering trans- 
form method cannot be applied to the corresponding boundary value prob- 
lem. The method requires very essential modification. Various ideas to 
extend the ISM to the initial boundary value problems are suggested in [12], 

HSMEMEI- 

In [5], [6] an effective tool to search integrable boundary conditions has 
been proposed based on some special involutions of the auxiliary linear prob- 
lem. This method (below for the sake of convenience we refer it as method 
of involutions) can be applied to integrable equations in both 1+1 and 1+2 
-dimensional cases. In this article we show that the method of involutions 
allows one to extend the concept of integrability to boundary value problems 
on bounded regions enclosed by any closed smooth curve. 

Let us explain briefly the approach we use. We call boundary value prob- 
lem integrable if it admits a Lax pair. Because of this reason we look for 
boundary condition simultaneously with its Lax representation. The start- 
ing point is to make a correct assumption about the possible form of the 
Lax pair of the boundary value problem. Actually this Lax pair is made 
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up from several different Lax pairs of the original equation itself by gluing 
the eigenfunctions along the boundary by properly chosen additional bound- 
ary conditions. As examples we take Liouville equation and two-dimensional 
Toda lattice equation. To generate new Lax pairs we use point symmetries 
(involutions) which leave invariant the nonlinear equation under considera- 
tion but change its Lax pair. 

In the next section, as a trial example, we consider the Liouville equa- 
tion. We remind the definition of integrable boundary conditions and find 
an example of integrable boundary conditions on a circle. It is shown that 
the following nonhomogeneous Neumann problem 

1 1 

u rr + -u r + —u 9 9 = 8e u , (1) 
2 

u r \ r=a 

a 

is integrable, i.e. it admits a Lax pair (see the list at the end of the second 
section). 

In the third section we study the two-dimensional Toda lattice equation 
on a circular cylinder: r < a, < 6 < 2n, — oo < n < oo. Several types of 
integrable boundary value problems for this lattice are found by using the 
method of involutions. Let uo(n) = exp{-u(n) — u(n + 1)}. It is shown that 
the following boundary values problems are integrable: 

U rr + ~^U r + ^U de = Uj(n — l) — Uj(n), (2) 

u{n)\r =a = UnQ + g{9) + k(n), 




U rr + ~Ur + ^ U 06 = ^(n - 1) - Uj(n), (4) 
l 

(u r (n) + -ue{n))\r =a = g(0), 

(X 

It is remarkable that the boundary conditions contain arbitrary functions g{9) 
and k{n). The Lax pairs for the above integrable boundary value problems 
are given in the list at the end of the third section. 
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In the fourth section we consider periodicity closure constraints reducing 
Toda lattice to the sinh-Gordon and Tcsitcseika equations. It is shown that 
the following boundary values problems are integrable: 
for the sinh-Gordon equation 

1 1 

Prr + -p r + ^Pee = 4 sinhp, (5) 

(p r + -pe)\r=a = 0; 

a 

for the Tcsitcseika equation 

q rr + -q r + \qee = e 2q - e~ q (6) 
(q r + -qe)\r=a = 0. 

Ci 

In the fifth section we give a class of exact solutions of the Toda lattice 
on a circle with nonhomogeneous Neumann type boundary condition on a 
circle. 



2 Liouville equation 

In this section we concentrate on boundary value problems for elliptic equa- 
tions. Suppose that the boundary T of a domain D is parameterized by the 
equation x' = f(t') that introduces a local system of coordinates by taking 
t-axis along the tangent direction and x-axis along normal to the curve T. 
Suppose that the differential equation under consideration 

E{u) = (7) 

admits two different Lax representations. For the sake of simplicity we take 
them rewritten in terms of the new coordinates 

Y x = U(X,u,u x ,...)Y(X) 

Y t = V(X,u,u x ,...)Y(X) {6) 

and 

Y x = U(X,u,u x ,...)YCX) () 
Y t = V(X,u,u x ,...)Y(X). 1 ] 



4 



where A, A are spectral parameters. Now the equation of the boundary is 
of the form x = 0. We are looking for conditions that allow to relate the 
equations for t evolution along the boundary, since x is fixed. More precisely, 
we have the following definition 



Definition 1 A boundary condition 







(10) 



is integrable if there exists a matrix F(X, t, u, ...) and function h(X) such that 
on the boundary x = the function Y = F(X,t,u, ...)Y(X) is a solution of 
the equation Y t = VY for any solution Y of the equation Y t = VY with 
X = h(X), provided the boundary condition holds. 

If a boundary condition is integrable in the sense of the definition above 
this means that the corresponding boundary value problem admits the Lax 
representation consisting of the two Lax pairs ([H]) and defined on the 
domain D such that the eigenf unctions Y and Y satisfy along the boundary 
an additional boundary condition (Y — FY)\r = 0. 

To consider a circle as a boundary we use polar coordinates (r,8). So, 
the boundary is r = a. In polar coordinates the Liouville equation is 



1 1 

u rr + -u r + -^u ee 



8e u . 



It admits the Lax pair 



where x = r, t 



Y r = LY, Y e 
-■ 9, U = L, V = A and 



2A 

o u+i0 i 



+ Xe 



-ie 



AY, 



(11) 

(12) 



2A 



l Ur + Tr Ue 



2X 



4 4r 



2A 



- Ae 



-w 



(13) 



i+i.e 



A 



IT 



2A 



- Ae 



-id 



u+i6 



2A 



1 i \ 

~4 Ur ~^ Ue 



V 2A 



2A 



+ Ae" 



(14) 
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To obtain a second Lax representation we use the Kelvin transformation. 
The equation (TTTj) is invariant under the Kelvin transformation 



a - A 1 a 

r = — u = u + 4 In — . 

Under such transformation the Lax pair (1121) takes form 

Y r = LY, Y e = AY, 

where 



L 



\ 2a?X 
( 



2a 4 A 



2 

r u,. 



Xe 



-i6 



2a 4 A 



r u r 



ir 



4a 2 a 2 4a 2 



Ug 



r ir 
Ao? a? Aa? U ° 



— la 
A = — 
r 



r 4 e u+i6 

2a 4 A 



-Ae 



2a 4 A 



r A e u+i6 r 2 



Ae 



-id 



(15) 

(16) 
\ 

/ 



2a 4 A 



+ 



4a 2 + a 2 



17) 



4a 2 



r 4 e u+ld r 2 u r 



r iriig 



\ 2a 4 A 4a 2 a 2 4a 2 



2o 4 A 



Ae" 



/ 



We will look for a boundary condition under which there exist a trans- 
formation A = h(X) and non degenerate matrix F(X,9,u) such that Y(X) = 
FY(h(X)) will solve the Y e (X) = A\ r=a Y(X) for every solution Y(X) of the 
equation Y e (X) = A\ r=a Y(X). 

Lemma 1 The integrable boundary condition is given by 

-2 



and there are two choices for the matrix F and the function h 
(i) 

h(X) = A, F= ' 1 ° 



(ii) 



h(X) = -A, 



F 



1 
- 



(19) 



(20) 



(21) 
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Proof. Let Y(X) satisfies the equation Yg = AY. On the boundary r = a, 
Y = FY(h(X)) has to satisfy Y e = AY. Substituting Y = FY(h(X)) into 
Yg = AY and using Yg = AY for Yg(h(X) we obtain 

±F - A(X)F + FA(h(X))\ f (MA)) = (22) 
The above equality holds if 

±F = A(X)F - FA(h(X)) (23) 

We have an equation for the unknown matrix F and function h(X). To 
solve the boundary condition ffTUl) with respect to u r we let u r = G(9,u,Ug). 
Assuming that F does not depend on ug and differentiating ( 1231) twice with 

respect to ug we obtain „ = 0. That is, u r = gi(u,9)ug + g 2 {u,9). We 

du z e 

substite the above expression for u r into the (1231) and let 



/ Mu,X,0) f 12 (u,X,9) \ 

V M(u,x,e) f22(u,x,0) J ' 1 } 

Separating terms with ug and without ug in (|23|) we obtain two sets of equa- 
tions. We write the first set of equations, terms with ug, as 

lf = Pf, (25) 
where / = (f n , f 12 , / 2 i, / 2 2) T and P is a matrix 



'/.a 



(26) 



-(f-i) ° ° (-t-a 

V o (-f-i) -(t-^) o y 

We write the second set of equations, terms without Ug, as 

£/ = «/, (27) 



where Q is a matrix 



r e u + i0\ 
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2A 



V2 

1 



V 



e"+' 9 _ g2 
2A 4 



£2 _ e u + iB 
' 4 2A 



-/i — V 



2h(X) J 



92 
4 



2A 



2h(X) 



-jJL + V 



(2i 



/ 



with 



2A - 2MA)- A ( A »^ m " 
The equations (125]) and (1271) must be compatible. This leads to the fol- 



92 1 
4 a 



lowing compatibility condition 



(P*-Q« + [P,Q])/ = 0, (29) 

where [P, Q\ is a commutator of P and Q. The matrix (Pg — Q u + [P, Q]) is 
nonzero. To have nonzero solution / the determinant of (Pg — Q u + [P, Q]) 
must be zero. It gives the following equality 

(h 2 (X)(a 9l (u, 6) - if - (a 9l (u, 9) + z) 2 A 2 ) 2 = 0. (30) 

The above equality holds if either 

(1) h(X) = p\ and g x = |gg where (3 G R\{-1, 1} or 

(2) h(X) = A and g x = or 

(3) h(X) = -X and g x = 0. 

One can show that in the case (1) there is no vector / to satisfy equations 
(1251) and ([27}. In the case (2) one has the only solution / = q(X) (1,0,0, 1) T 
if g 2 = =1, That gives the boundary condition ([191) . function h and matrix 
F given by (T2D1). 

The case (3) is similar to the case (2) and gives same boundary condition 
(jl"9l . function /i and matrix F given by (I21"I) .D 



From the above lemma we have the following integrable boundary value prob- 
lem with corresponding Lax pairs (we have two Lax pairs for the problem). 
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1 1 

r < a u rr H — u r H — -uee = oe , 

F r (A) = L(A)V(A), Y r (X) = L(X)Y(X), 

Y e (X) = A(X)Y(X), ana Y e (X) = A(X)Y(X), 
2 

I Qi IXf y — d 

a 

Y = FY 

where F is given by (120]) . L is given by (TT3"j) . A is given by (THI) and L 
is given by (TT7|) . A is given by fTTHT) . 

1 1 

r < a u rr + -u r + —Uee = 8e , 

Y r (X) = L(X)Y(X), Yr (-X) = L(-X)Y(-X), 

Y e (X) = A(X)Y(X), and Y e (-X) = A(-X)Y(-X), 
2 

T (X Ur\r= a 

Y{X) = FY(-X) 

where F is given by fl2T|) . L is given by f|T3|) . A is given by ffT^l) and L 
is given by f|T7|) . A is given by f|T8l) . 



3 Two-dimensional Toda Lattice 

We make the same assumption, as in the case of the Liouville equation, for 
the coordinates. Hence, boundary is given by x = 0. Again we suppose 
that the differential equation under consideration admits two different Lax 
representations 

Y X = UY . Y x = UY , Q1 , 

Y t = VY and £=VY (31) 



For the two dimensional Toda lattice equation U, V, U and V in (1311) are 
linear operators. 

Definition 2 A boundary condition 

tt(u) = (32) 

is integrable if there exists a linear differential operator A such that on the 
boundary x = we have Y = AY is a solution of Y t = VY for any solution 
YofYt — VY, provided the boundary condition holds. 
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To consider a circle as a boundary we use polar coordinates (r,8). So, 
the boundary is r = a. The two dimensional Toda lattice equation in polar 
coordinates becomes 

1 1 

u rT + -u r + -^u ee = uo{n - 1) - u(n), (33) 
where w(n) = exp(u(n) — u(n + 1)). The above equation admits a Lax pair 
e ie If i \ 



e -iB 



u(n- l)0i(n- 1), (34) 



ViA n ) = — + 1) - 7T u r(n) Ufl(n) ¥?i(n)+ 



r 

— is 

+!^_ w (n- l)^(n-l) (35) 

To obtain other Lax representations we use symmetries of the equation fl33l) . 

1. Reflection on # 

= -5; (36) 

2. The Kelvin transformation 

f = — , u = u + Anln (37) 

3. Reflection on n 

u = -u(-n). (38) 
Using the transformation ( 1361) we obtain the following Lax representation 

¥2,r(n) = —Lp 2 (n + 1) - - lw r (n) + -u^(n) j <£ 2 (n)- 

-—ou(n-l)ip 2 (n-l), (39) 
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¥2,e{n) = — - — ip 2 (n + 1) - — ( u r {n) + -ug(n) J (p 2 (n) + 

id 

+ l l?-u{n-l)ip 2 {n-l) (40) 

Using the Kelvin transformation (1371) we obtain the following Lax represen- 
tation 

g /* f 2,f "\ 

V3,r(«) = -^-^(n + 1) - - I — ttr(n) + 4n— ^ 4 (n)- 

Z V CJ Ci Ci / 



2a 4 



w(n-l)p 3 (n-l), (41) 



¥>3,*(n) = -^r<P3(n + 1) - — ( — Ur(n) + 4n— - ) <^ 3 (n)+ 



2r 2r V a 2 a 2 a 2 

w(n-l)^ 3 (n-l). (42) 



• 3 — ifl 

ire 



2a 2 

Using the transformations ( 1381) we obtain the following Lax representation 



e ie If i 

<P4,r( n ) = -TT^in ~ 1) - - -u r (n) + -u (n) ) <p±(n)- 



2 r y ' 2 V v ' r 



e -iB 



u(n)<p 4 (n + l), (43) 



<fi4:,e( n ) = -^-^i 71 - 1) _ ~2 ( - u r(n) + j Lp 4 (n) + 

—id 

+ !!^_ w ( n )^ 4 ( n+ 1). (44) 

According to Definition 2, to obtain the integrable boundary conditions we 
relate the equations for 9 evolution of the above Lax representations, on the 
boundary r = a. We consider the case when the eigenfunctions are related 
by multiplication operator ipt = A(9, n,u, . . . ) • (fj. 

It turns out (see the lemma [7J that that Lax pairs corresponding to the 
Kelvin transformation fl37j) and the symmetry (138]) are gauge equivalent. 
A solution of (1421) transforms to a solution of (|44p without any boundary 
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conditions. So, for boundary value problems ([2D and (J3J) we have two possible 
Lax pairs. 

In the lemma [2] we derive the Lax pair for the boundary value problem (T5|). 

Lemma 2 Let (pi(n) be a solution of the equation ( TffSp then on the boundary 
r = a a function ^{ti) = A ■ (fi(n), where 

A = e 2m+g ^ e \ g(9) is an arbitrary function of 6, (45) 

is a solution of the equation [JW provided the following boundary conditions 

u(n) = 2%n9 + g(9) + k(n), k(n) is an arbitrary function ofn, (46) 

holds (for all n). 

Proof. On the boundary r = a we substitute ^>2{ n ) — A(n, 9,u, . . .) ■ <pi{n) 
into the equation (HOI) and use (1551) for (^^(n). The resulting equation holds 
if the coefficients of <pi(n+ 1), fi{n) and tpi(n — 1) are zero. Thus we obtain 
the equations 

%ae %e . , . iae~ l " 



2 



A(n) = — — A(n + 1) (47) 



, , , la I , , lUfj(n) \ . la , s iua[n \ ., . . , . 
4>M - y U-M - J = -y kW + -^J A(n) (48) 

^— w(n - l)A(n) = — w(n - l)A(n - 1) (49) 

From the equations (T4T1) . (|49p we have that A(n) = e 2l9 A(n — 1). Hence, 
A(n) = e 2l6n b(9) where b(9) is a function of # only. Substituting A(n) = 
e 2l9n b(9) into the equation (|48p we obtain 

b e + (2m - u 6 {n))b = 0. (50) 

Since the function b does not depend on n we have that the coefficient of b in 
the above equation does not depend on n, so u$(n) = 2in + h{9). Integrating 
with respect to 9 we obtain the boundary condition (j4"6"l) . where k is an 
arbitrary function of n and g{9) = J h(9)d9. Then solving the equation 
(|48l) . assuming that the found boundary condition holds, we obtain A = 

e 2in+fg(9)d9^ the express i on ggj) f or A Q 
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In a similar way, from the next lemmas we have Lax pairs for the boundary 
value problems fl3]) and (j3J. 

In the lemma [3] and the lemma H] we derive the Lax pair for the boundary 
value problem [3j 

Lemma 3 Let <pi(ri) be a solution of the equation If35*\) then on the boundary 
r = a a function ips{n) — A • <pi(n), where 

A = e m f 9 ^ de , g{9) is an arbitrary function of 9, (51) 

is a solution of the equation ffllfy provided the following boundary conditions 



In 

u r (n) = — + g(9) (52) 



holds. 



Proof On the boundary r = a we substitute (ps(n) = A(n) ■ <p\{ri) into the 
equation (1421) and use (135]) for (px^{n). The resulting equation holds if the 
coefficients of <fi(n + 1), tpi(n) and (pi(n — 1) are zero. Thus we obtain the 
equations 



iae 10 , , , i(i< 



} A(n) = —A(n+1) (53) 

. . . ia ( . . iu s (n)\ ... ia ( , s An iu e (n)\ 

M n ) - ~77 M n ) — Mn) = — u r {n) + — ^ A(n) 

2 \ a J 2 \ a a J 

(54) 

—id ■ —id 

w(n - l)A{n) = uj{n - l)A{n - 1) (55) 

From the equations ( 1531) . ( 1551 we have that ^4 does not depend on n. Hence, 
the coefficient of A in ( 1541 must be a function of only. This gives the 
boundary condition (|52|) . Then solving the equation (1541 . assuming that 
(|52|) holds, we obtain the expression (l5Tl) for A. □ 

Lemma 4 Let (fii(n) be a solution of the equation l[35\) then on the boundary 
r = a a function (pi(n) = A ■ (pi(n), where 

A = e 2i9n+u(n)+iafg(6)d9^ g ^ w m arbUmry f unctlon f Q ; ( 56 ) 

is a solution of the equation provided the following boundary conditions 

2n 

u r {n) = — + g{9) (57) 
a 

holds. 
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Proof. On the boundary r = a we substitute (p±{n) = 9,u, . . .) • (fi(n) 
into the equation (jUJ) and use (jl2"l) for (^^(n). The resulting equation holds 
if the coefficients of (p% (n + 1), (f2{n) and — 1) are zero. Thus we obtain 
the equations 

id ■ —id 

A{n) = % ^^u(n)A(n + 1) (58) 



2 y ' 2 

M n) - f (u,.(n) - ^) A(n) = -| 4- *4=>) A(„) (59) 

^|_ £4 ;(7i-l)A(n) = ^-A(n-l) (60) 

From the equations (|5^|) . (I6"U1) we have that A(n) = e~ 2ie ti;(n)yl.(n+l). Hence, 
A(n) = e 2i6n + u ( n )l)^ where F{9) is a function of 6* only. Substituting A{n) = 
e~ 2ten b(9) into the equation (1591) we obtain 

be — (2m — iau r {n))b = 0. (61) 

Since the function b does not depend on n we have that the coefficient of 
b in the above equation does not depend on n. This gives the boundary 
condition (1571) . Then solving the equation fl59l) . assuming that the found 
boundary condition holds, we obtain the expression f[56"j) for A. □ 

In the lemma [5] and the lemma [6] we derive the Lax pair for the boundary 
value problem HI. 

Lemma 5 Let (p2(n) be a solution of the equation HW then on the boundary 
r = a a function y?3(n) = A ■ ^>2{n), where 

e -2ien+iaf g (e)d8^ g(#) i s a n arbitrary function of 6, (62) 
is a solution of the equation fiJJfy provided the following boundary conditions 



u r (n) = — u e {n) + g{6) (63) 
a 



holds. 



Proof. On the boundary r = a we substitute <^>z{ n ) — A(n, 0, u, . . . ) ■ <f2(n) 
into the equation (l4"2l and use (14"0"1) for ip 2 ,e(n). The resulting equation holds 
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if the coefficients of p%{yi + 1), p>2{ n ) and (p2{n — 1) are zero. Thus we obtain 
the equations 



iae 10 . , , toe 



2 



A(n) = A(n + 1) (64) 



a i \ i a ( i \ iu B (n)\ ,, . m / , . 4n iu e (n)\ , 
2 V a / 2 V r a 



if) ■ —if) 

iae „ , , zae 



(65) 



2 A(n) = —j-A(n-l) (66) 

From the equations fl64l) . fl66l) we have that A(n) = e _2i6 'y4(n — 1). Hence, 
A(n) = e~ 2l6n b{6) where b(8) is a function of 6 only. Substituting A(n) = 
e~ 2l6n b{6) into the equation ( 1651) we obtain 

be — ia(u r (n) H — Ue(n))b = 0. (67) 

Since the function 6 does not depend on n we have that the coefficient of 6 in 
the above equation does not depend on n. This give the boundary condition 
(|63|) . Then solving the equation (|65|) . assuming that the found boundary 
condition holds, we obtain the expression (1621) for A. □ 



Lemma 6 Let </2 2 (n) &e a solution of the equation Hlty then on the boundary 
r = a a function tp±(n) = A ■ P2{ n ), where 

A = e «(»)+/ g (p) ls an arbitrary function of 6, (68) 
is a solution of the equation U^f} provided the following boundary conditions 

u r {n) = -\ e {n) + g{9) (69) 

holds. 

Proof. On the boundary r = a we substitute tp^n) = A(n, 9, u, . . . ) • ^2(71) 
into equation (jHJ) and use (j4"01) for (p 2 ,e( n )- The resulting equation holds if 
the coefficients of + 1) 3 ^(n) and ^2(^1 — 1) are zero. Thus we obtain 
the equations 

—if) ■ —if) 

iae , , , iae 



A(n) = — — w(ra)j4(ra + 1) (70) 
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Mn) - | U(n) + ^ A(n) = -f f -„ r (n) + ^ ] -4(„) (71) 



iO ■ i0 

iae , , . , , zae 



2 w(n-l)i4(n) = — A(n-1) (72) 

From the equations (1701) . (1721) we have that A(n) = u(n)A(n + 1). Hence, 
A(n) = e u{n) F{8) where F{8) is a function of 8 only. Substituting A(n) = 
e u(j ^b{8) into the equation flTTl) we obtain 

+ (—iau r (n) + ug(n))b = 0. (73) 

Since the function b does not depend on n we have that the coefficient of b 
in the above equation does not depend on n. We obtain the boundary condi- 
tion (|69|) . Solving the equation (ITT]) and assuming that the found boundary 
condition holds, we obtain the expression (1681) for A. □ 

In Lemma [7J we show that the Lax representations corresponding to the 
Kelvin transformation (l3Tj) and the symmetry (1381) are equivalent. 

Lemma 7 Lei (^(n) 6e a solution of the equation ( f^I| ) i/ien on t/ie boundary 
r = a a function (pi(n) = A ■ (p 2 (n) , where 

is a solution of the equation fi44\ )- 

Proof. On the boundary r = a we substitute ^(n) = A(n, 8, u, . . . ) • ^(ft-) 
into the equation (jUJ) and use (HT)1) for ^^(n). The resulting equation holds 
if the coefficients of Lp^(n + 1), 923(71) and 933(7?, — 1) are zero. Thus we obtain 
the equations 

^(n) = ^( n )i( n + l) (75) 

Mn) ~ * (-u r (n) + ^ - ™°&) A(n) = f-^n) + ^) A(n) 

(76) 



2 \ v ' a a / v ' 2 



— iO ■ it) 

iae , \ * / \ iae 



y u(n-l)A(n) = —A{n-l) (77) 

From the equations (!70|) . (1721) we have that A(n) = e~ 2lk u(n) A(n+1) . Hence, 
A(n) = e 2ien + u( , n )l){Q^ where b(8) is a function of 8 only. Substituting A(n) = 
e 2ien+u(n) b ^ intQ the equat j ori (jjjj) we obtain 

& e = 0. (78) 
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Hence the function b is a constant. This gives us the expression (I68p for A. 
□ 



From the above lemmas we have the following list of integrable boundary 
value problems with corresponding Lax pairs. Some of the integrable bound- 
ary value problems admit two different Lax pairs. We give both Lax pairs in 
the list. 

The list of Lax pairs for two dimensional Toda lattice 

• r < a u rr H — u r H — -u ee = ou(n — 1) — u)(n), 

(fi, r {n) = E/i<pi(n), <P2,r( n ) = U2<f2{n), 

<Pi,e(n) = Vi<fi(n), P2,e( n ) = ^^(n), 

r = a u(n) = 2ind + g{9) + k(n), 

where action of operator U\ is given by fl34l) . V\ is given by fl35l) and 
U2 is given by ( |39l) . V2 is given by rtlDjl . 

• r < a u„. H — u r H — -m^ = cj(n — 1) — u;(?t,) 

^i, r (n) = [/i</?i(n), ^ ^3,r(«) = U 3 tp 3 (n), 
<Pi,e( n ) = Viipi(n), f3,e(n) = V 3 (p 3 (n), 

r = a u r [n) = h g{9), 

a 

ip 3 = e^IaWde^ 

where action of operator U\ is given by (Ell), V\ is given by ( 1351) and 
U 3 is given by (|4T|) . V3 is given by (Pf2j) . 

1 1 , _ , . 

• r < a u rr H — u r H — jrw^ = uj(n — 1) — u;(n) 

^i, r (n) = U x y x {n), ^ ^ 4 ,r(«) = U 4 cp 4 (ri), 
<Pi,e{ n ) = Vi^iO 1 )) ^4,0 (n) = V^O), 

r = a it r (n) = h 

a 

yj 4 = e 2ien+u(n)+iaJ g{8)dd ( ^ i 

where action of operator is given by (1311) . Vi is given by (1351) and 
U4 is given by (|43|) . V4 is given by (jHJ). 

1 1 , . . . 

• r < a u rr H — u r H — riise = cum — 1) — wm) 
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¥V(n) = U 2 <p 2 {ri), &nd wfa) = U 3 <p 3 (n), 

(P2,e(n) = V 2 (p 2 (n), <P3,e( n ) = V 3<Ps(n), 

i 

r = a u r (n) = — ug(n) + g(9), 

a 

tp 3 = e -2i8n+ia J g(9)de ^ 

where action of operator U 2 is given by fl39l) . V 2 is given by fHUl) and 
[7 3 is given by ffldj) . V 3 is given by f|42|) . 

1 1 

• r < a u rr H — u r H — = ^(^ — 1) — uo(n) 

(p 2<r (n) = U 2 ip 2 (n), tpi, r {n) = U±tp±{n), 

i 

r = a u r (n) = — ug(n) + g(9), 

where action of operator U 2 is given by (1591 . V 2 is given by (HOI and 
U4 is given by (H3|) . V4 is given by (T44|) . 



4 Reductions of two dimensional Toda lattice 
equation 

In this section we obtain integrable boundary conditions for the sinh-Gordon 
and Tcsitcseika equations as reductions of integrable boundary conditions of 
the two dimensional Toda lattice equation. 

To reduce the two dimensional Toda lattice equation to the sinh-Gordon 
equation we put periodicity condition u{n) = u{n + 2) for all n, where u 
satisfies ( 1551) . Then for p = u(0) — u(l) we have 

1 1 

Prr + ~Pr + ~^Pee = 4 sinhp, (79) 

the sinh-Gordon equation in the polar coordinates. Only the boundary con- 
dition of the problem (T4J) is consistent with periodicity constraint u(n + 2) = 
u{n). It gives 

p r + - Pd = (80) 
a 

on the boundary r = a. Evidently by changing p = iv we get v rr + -v r + 
\v ee = 4sint> and v r + -v e = 0. 
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To reduce the two dimensional Toda lattice equation to the Tcsitcseika 
equation we put u(n) = u(n + 3) and u(n) = —u(2 — n). Then for q = u(0) 
we have 

q„ + -q r + —q ee = e 2q - e~ q (81) 

the Tcsitcseika equation in polar coordinates. Again only the boundary con- 
dition of the problem (jj]) is consistent with periodicity constrant u(n) = 
u(n + 3) and u(n) = —u(2 — n) . It gives 

q r + -q e = (82) 
a 

on the boundary r = a. 



5 Some solutions of the boundary value prob- 
lem 

In this section we give an example of solutions for the special case of the 
boundary value problem 

11 2n 

u rr + -^u ee + -u r = to(n - 1) - u(n), u r (n)\ r=a = h g(9) (83) 

r z r a 

where u(n) = exp(u(n) — u(n + 1)). We assume that g{6) = and look 
for spherically symmetric solution. That is u is a function of r only. The 
boundary value problem ( 1831) reduces to 

1 2n 
u rr H — u r = u(n — 1) — u(n), u r (n)\ r=a = — . (84) 
r a 

Let us introduce new variables t = lnr and v(n,t) = u(n,r) — 2nlnr. Then 
the boundary value problem (IMj) becomes 

v tt = u>(n-l)-u(n), u t (n)\ t=0 = 0, (85) 

where u>(n) = exp{v(n) — v(n+l)}. As solutions of the above boundary value 
problem we can take even solitons of Toda lattice equation in one dimension. 
Following [18], (see pp. 494-498), the general iV-soliton solution is given in 
terms of the data {c, Zj, 7j} such that 
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I. The quantities Zj lie in the interval — 1 < Zj < 1 and are pairwise 
disjoint. 

ii. e - c =nf=i4 

TV 

III. The quantities m,-(0) = , 3 . , where a(z) = TTsignz,- — and dot 

means derivative with respect to z, are positive. 

The iV-soliton solution is given by 

, N , det M(n,t) 

"'"■"-^'■' detM^-M) (86 > 

. , . . Jm i (t)m j (t)(z i z 1 ) n+l 

where M(n, t) is a matrix with entries Mij(n, t) = 5ij + 



and mj(t) = — — ^ — -, i, j — 1, . . . N. 



1 Zi Zj 



a[Zj i 



The even solitons are describe by the following lemma. 

Lemma 8 Let N = 2k and the data Zj, ^j, j — 1 . . . N satisfy z^ = —z^-i+x, 
7i = — 7/v-i+i, i — 1 . . . k. Then the N soliton solution is even function 
oft. 

Proof. With our choice of the initial data the elements of matrix M(n,t), 
that are symmetric with respect to the "center" of the matrix, are equal. If 
t is changed to —t then every element of M(n, t) is replaced by the element 
symmetric to it with respect the "center" of the matrix. Hence determinant 
of M(n,t) is equal to the determinant of M(n,—t) and v(n,t) = v(n,—t). 
From v(n,t) = v(n, —t) it follows that v'(t)\ t= o = O.D 

We give an example of the solutions described in the above lemma. For 
N = 2 we put Z\ = z , z 2 = —z , c = — 4\nz and 71 = —70, 72 = 7o where 
< zq < 1 and 70 > 0. Then the data satisfy the conditions I, II, III and 
the conditions of the lemma. With such data one has the following solution 
of ([85D 

{ , ) 1 + 70(1 + ^)^- 1 cosh[(, - ^1+70^ ( j 
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Hence the boundary value problem (I84p has the following solution 



u(n, r) = c + In 




An+4 



An 




8. 



■8 



6 Conclusion 

In the present paper we apply the method of involutions to boundary value 
problems for soliton equations on bounded regions. As illustrative models we 
consider Neumann type boundary value problem on a circle for the Liouville 
equation and initial boundary value problem for the two dimensional Toda 
lattice equation. The Lax representations for the boundary value problems 
are represented. We considered some reductions of the integrable boundary 
value problems in the case of the two dimensional Toda lattice equation. We 
also constructed a class of solutions satisfying one of the found boundary 
conditions. 
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